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Abstract 

The IceCube Neutrino Observatory is a kilometer-scale detector 
currently under construction at the South Pole. The full detector will 
comprise 5,160 photomultipliers (PMTs) deployed on 86 strings from 
1.45-2.45 km deep within the ice. As of the austral summer of 2009-10, 
73 out of the total number strings have been deployed, and the detector 
is reaching its final construction phase. A dense sub-array of 6 strings 
in the center of the detector (DeepCore) has been already installed for 
enhancing the sensitivity to low energy neutrinos. The IceCube de- 
tection principle is based on the measurement of the Cherenkov light 
induced by ultra-relativistic muons and showers produced by neutrino 
interactions in the target matter of the detector. The main scientific 
goal of the IceCube experiment is the detection of astrophysical neu- 
trinos that will help to understand and settle the unresolved questions 
about the origin and nature of cosmic rays. In this contribution we will 
present the latest results of the experiment concerning the search for 
neutrino point sources using the experimental data taken during 2008- 
09 where the detector was operated with a 40-string configuration. The 
results of the analysis for steady individual neutrino sources as well as 
the stacking analysis from different catalogs will be presented. 



1 Introduction 

The IceCube Neutrino Observatory is a neutrino telescope designed to de- 
tect high energy astrophysical neutrinos with energies > 100 GeV. Such 
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an observation could reveal the origins of cosmic rays (CR) and the possi- 
ble connection to shock acceleration in Supernova Remnants (SNR), Active 
Galactic Nuclei (AGN) or Gamma Ray Bursts (GRB). An advantage of neu- 
trino astronomy over other experimental observations such as gamma ray 
astronomy is the possibility to offer insight about the most energetic and 
dense parts of the Universe since neutrinos are weakly interactive particles 
that can travel through matter without interacting. Nonetheless, this ad- 
vantage also makes the neutrino detection a technological challenge since a 
large amount of volume is required to observe neutrino interactions. The 
IceCube Neutrino Observatory uses the Antarctic ice as a detection volume 
where a neutrino interaction will produce a muon (or other charged lepton) 
as a result of the charged current interaction. Deep Antarctic ice is also 
transparent which allows the induced Cherenkov light by the muon to prop- 
agate. In order to detect this Cherenkov light, IceCube will instrument a 
cubic kilometer of clear Antarctic ice sheet underneath the geographic South 
Pole with an array of 5,160 Digital Optical Modules (DOMs) deployed on 86 
strings from 1.45-2.45 km deep. The DOMs are spherical, pressure resistant 
glass housing, each containing a 25 cm diameter Hamamatsu photomulti- 
plier tube and the electronics for waveform digitization. The IceCube con- 
figuration also includes a denser array with 6 strings equipped with higher 
quantum efficiency DOMs, the so called DeepCore detector, in order to in- 
crease the sensitivity to low energy neutrinos (< 100 GeV). A surface array 
for observing extensive air showers of cosmic rays is also being installed 
(IceTop). IceCube construction started with a first string installed in the 
2005-6 season (Achterberg et al. 2006) and will be completed in the austral 
summer of 2010-11. 

In this paper we describe the point source analysis of the data corre- 
sponding to the year 2008-9 when IceCube configuration consisted in 40 
deployed strings. The IceCube 40-string configuration performance as well 
as the event selection criteria followed in the analysis are described in sec- 
tion [2 Section [3] explains the methodology used for the search of point 
neutrino sources and stacking analysis. The search strategies applied to the 
IceCube data and their results are discussed in section [H Conclusions are 
shown in section [5j 

2 Event selection and detector performance 

As mentioned before, forty strings of IceCube were operational from April 
2008 to May 2009 with ~ 90% duty cycle after a run selection based on the 



stability of the detector. The number of trigger events of the IceCube de- 
tector for the 40-string configuration livetime was of the order of ~ 3 x 10 10 
events. Among these events, 99.99999% are muons produced by the impact 
of primary cosmic rays with the atmospheric nuclei. This number of trig- 
gered events can be reduced to ~ 8 x 10 8 events using an online filtering sys- 
tem based on the quality of the track reconstructions and energy estimators. 
The total livetime of the 40-string configuration was 375.5 days, leading to 
final atmospheric neutrino rate of 40 atmospheric neutrinos per day, which 
comprises the irreducible background of a point source analysis. The de- 
tection principle of neutrino telescopes normally uses the Earth as a shield 
for up- going atmospheric muon tracks. By selecting only reconstructed up- 
going tracks we can reject part of these atmospheric muons. However due 
to the huge amount of down-going atmospheric muons reaching the detec- 
tor, even a small fraction of mis-reconstructed events can contaminate the 
up-going sample. For that reason tight cuts in the quality of the track re- 
construction and track-like parameters such as the reduced likelihood of the 
track fit and the directional width of the likelihood space around the best 
track fit are applied. Also a cut in the likelihood ratio between the best up- 
going and down-going track solution is applied together with a requirement 
that the event's set of hits can be split into two parts which both reconstruct 
as nearly up-going. These cuts in the up-going region (northern sky from 
IceCube location) define a final data sample with sensitivity to point sources 
optimized for sources of neutrinos in the TeV-PeV energy range. 

With the 40-string configuration we also extended the point source anal- 
ysis to directions above the horizon in the down-going region of the sky 
(southern sky from the IceCube location). This extension of the point 
source analysis is possible by rejecting five orders of magnitude of atmo- 
spheric muon background events by applying cuts in the energy estimator 
in order to discriminate the large amount of these atmospheric muons from 
a possible neutrino signal with a harder spectrum. After track-quality se- 
lections, similar but tighter than for the up-going sample, a cut based on 
an energy estimator is made until a fixed number of events per steradians 
is achieved. Because only the highest energy events pass the selection, sen- 
sitivity in part of the sky is primarily to neutrino sources at PeV energies 
and above. 

Both the angular resolution and effective area of the detector depend on 
the event selection. Since the two samples (up-going and down-going events) 
are defined by different quality cuts, the detector performance is different 
for both parts of the sky. Figure [1] shows the neutrino point spread function 
of the IceCube 40-string configuration for up-going events compared to the 
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Figure 1: Cumulative point spread function. The horizontal axis represents the 
angle between the true neutrino direction and the reconstructed muons. The ver- 
tical axis is the fraction of events with an angular distance less or equal to the 
corresponding value on the horizontal axis. 
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Figure 2: Angular resolution, defined as the median of the point spread function, 
for selected events in the up-going region (blue) and down- going region (red). The 
dashed area indicates a 10% area. 



past 22-string configuration and the final IceCube detector. The angular 
resolution, defined as the median of the point spread function, is shown in 



figure [2] for both the up-going and down-going region for the final IceCube 
40-string sample. After the event selection described in the previous section, 
the final data sample consists of 36900 events where approximately one third 
are up- going events and two thirds down- going events. 




Figure 3: Equatorial sky map of 36900 events in IceCube 40-strings after the anal- 
ysis cuts. 



3 Methodology 

The final sample of a point source analysis consists of a set of arrival direc- 
tions and energy estimators. Figure [3] shows how the final arrival directions 
are distributed on the whole sky. The aim of a point source analysis is to 
identify an anisotropy on these arrival directions by looking for an accu- 
mulation of events in a given direction which will indicate the presence of 
an astrophysical point source emitting neutrinos. In order to perform such 
an analysis in IceCube we use an unbinned likelihood method (Braun et al. 
2008). This method has substantial improvements with respect to the stan- 
dard binned analysis. The method accounts for individual reconstructed 
events uncertainties and energy estimators and compares the random dis- 
tributions of atmospheric neutrino background for the northern sky and 
atmospheric muon background for the southern sky to that produced by a 
signal coming from a point neutrino source. The method models the data 
as a two component mixture of signal and background. A maximum likeli- 



hood fit is applied to determine the relative contribution of each component. 
Given N events in the final sample the density distribution of the i th event 
is given by 



where Si is the probability density distribution for the signal hypothesis and 
Bi for background. The parameter n s is the number of signal events and is 
one of the free parameters of the likelihood maximization together with the 
spectral index, 7, of the signal spectrum distribution. 

For each tested direction in the sky, the best fit according to the like- 
lihood maximization is found yielding the best estimate of h s and 7. The 
logarithm of the likelihood ratio between the best-fit hypothesis, £(n s ,7), 
and the null hypothesis of only background, C(n s = 0) is the test statistic of 
the method. The significance of the result is evaluated by comparing the test 
statistic with a distribution obtained by performing the same analysis over 
a sample of only-background data sets. These only-background datasets are 
obtained by scrambling real data in right ascension. The Earth rotation 
guarantees the uniform exposure in right ascension and uniform background 
per declination band, hence we can randomize the events in right ascension 
while keeping all other event properties fixed. Nearly vertical events are 
excluded from the randomization since polar regions in the sky cannot be 
scrambled in right ascension. 

In all-sky analyses the whole sky is scanned in steps of 0.1° x 0.1° forming 
a grid with bin size much finer than the angular resolution. The hottest 
spot in the grid is considered to be the point with the highest significance. 
The final significance is determined by the fraction of scrambled data sets 
containing at least one grid point with a log likelihood ratio higher than 
the one observed in the data. This fraction is the post-trial p-value for the 
all-sky search. 

A significant penalty is paid in the all-sky search due to the number of 
effective trials due to scanning the entire sky. Therefore an additional search 
usually performed in point source analysis is to select a number of locations 
a priori based on candidate sources. The likelihood is maximized on each of 
these locations instead of over the entire sky reducing the number of trials. 
The most significant result from this candidate source list is also compared 
to the distribution of randomized scrambled data yielding the final post-trial 
p-value. 

A third method of stacking analysis is also presented in this paper where 
instead of considering individual sources the cumulative signal from a collec- 
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tion of similar sources is considered. The stacking analysis is a well-known 
technique in astronomy that has been already applied to neutrino astron- 
omy (Abbasi et al. 2009). Stacking multiple sources is an effective way to 
enhance the discovery potential of similar sources when the source fluxes 
are too faint to yield a discovery individually. By considering the cumula- 
tive signal for all sources emitting together, we can increase our discovery 
potential and set upper limits on interesting catalogs of neutrino source can- 
didates as a whole. The stacking analysis requires a small modification to 
the unbinned likelihood analysis. The signal density distribution, <Sj, from 
equation Q] now refers to the sum of individual sources. Assuming that M 
is the number of sources to stack, the signal density distribution can be 
rewritten as follows 

where Wj is a relative theoretical weight of the j th source in the catalog 
and Rjd) is the relative detector acceptance for a source with spectral 
index 7. As in the all-sky search and the fixed source search, the only two 
free parameters in the likelihood maximization are the signal contribution, 
n s , and the spectral index, 7. The weights, Wj are correlated with the 
theoretical expected neutrino luminosity from that given source based on 
gamma observations, a uniform model will assume all Wj = 1. The detector 
acceptance, Rji'y), is tabulated using simulation. 



4 Searches and results 

In the previous two sections we discussed the methodology followed in the 
point source analysis. With the method described we performed five different 
searches using the data of the IceCube 40-string configuration: 

• an all-sky search using the standard point source analysis. 

• a dedicated search based on a list of candidate sources to reduce the 
penalty from the number of trials. 

• a stacking search using the TeV source list from Milagro as a catalog. 

• a stacking search using a catalog of nearby clusters of galaxies. 

• a stacking search using a catalog of Stardust galaxies. 



4.1 All-sky search and candidate source list results 

The first search is a scan looking for the most significant point in the whole 
sky. As mentioned before, with the 40-string configuration the search was 
extended to the southern sky therefore covering declinations from —85° to 
85°. Figure H] shows the significance skymap for the all-sky search. The 
most significant deviation from the only background hypothesis is located 
at coordinates 113.75° r.a., 15.15° dec. The pre-trial estimated p-value of 
the maximum log likelihood ratio at this location is 5.2 x 10~ 6 equivalent to 
a single sided sigma of 4.5. 1817 out of 10,000 scrambled data sets have an 
equal or higher significance somewhere in the sky resulting in the post-trial 
p-value of 18% well compatible with a background fluctuation. 




Figure 4: Equatorial skymap of pre-trial significances (p-value) of the all-sky point 
source. The galactic plane is shown as the solid black curve. 

The candidate source list is formed by 39 a priori selected locations in 
the sky (see table [T]) . The most significant results comes from the location 
of PKS 1622-297 with a final post-trial p-value of 62% determined as the 
fraction of scrambled data sets with at least one source with an equal or 
higher significance. This result is also compatible with the only-background 
hypothesis. 
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Table 1: Results for the source candidate list. $go is the upper limit of 
the Feldman-Cousins 90% confidence interval for an E~ 2 flux, i.e.: d<fr/dE < 
$90 10 _12 TeV -1 cm~ 2 s -1 (i±;/TeV)~ 2 . n s is the best-fit number of signal events; 
when n s > the (pre-trial) p-value is also calculated. Ni° is the actual number of 
events observed in a bin of radius 1°. The background event density at the source 
declination is indicated by the mean number of background events B\° expected in 



4.2 Stacking search results 

For the stacking analysis three different catalogs were chosen. The first 
search is based on the observations of the Milagro collaboration which re- 
ported TeV gamma-ray emission from 16 sources after correlating it with the 
Fermi Bright Catalog list (Abdo et al. 2009). These sources are promising 
candidates for neutrino emission (see Beacom & Kistler 2007 as an exam- 
ple). An additional source was confirmed to have TeV emission, the MGRO 
J1852+01. We perform a search based on 17 sources observed in TeV gamma 
rays by Milagro using an equal theoretical weighting from equation [2j The 
second stacked search is performed using 127 starbust galaxies. Starbust 
galaxies have a dense interstellar medium (ISM) and high star formation 
rates. The elevated CR rate due to the high number of Super Nova Rem- 
nants (SNR) and the dense ISM leads to a potential neutrino emission from 
these starbust galaxies (see Becker et al. 2009). 

Cluster of galaxies are also a potential source of high energy protons 
and therefore neutrinos. The third stacked analysis uses 5 clusters of galax- 
ies using the model predictions in Murase et al. 2008 as weights in the 
likelihood. 
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Table 2: Results for the stacked source searches. The p- value indicates the proba- 
bility of background to produce the same or higher observed significance. The flux 
is the upper limit (90% C.L.) set to the corresponding model. 



The results for the stacking analysis are summarized in table [2j No 
significant excess has been found and therefore upper limits for a given 
model have been calculated. In the case of Milagro sources and cluster of 
Galaxies, the upper limit shown on table[2]are calculated assuming a uniform 
theoretical weighting and an E~ 2 neutrino emission. The upper limits for 
the starbust galaxies stacking analysis is shown for a model with weighting 
proportional to the far infrared (FIR) flux at 60 [im since the FIR and radio 
emission from these sources are correlated with hot dust and a high star 
formation rate given in Table A.l in Becker et al. 2009. 



5 Conclusions 

The IceCube observatory has analyzed the data corresponding to the year 
2008-9 when the detector configuration consisted of 40 strings. The analysis 
of time integrated point sources showed no evidence of a neutrino signal 
and upper limits have been calculated. For an E~ 2 signal of Vn the cal- 
culated upper limits from this work are shown in figure El The predicted 
sensitivities of IceCube and the ANTARES neutrino telescope located in the 
Mediterranean Sea are also shown. The sensitivity is to be interpreted as the 
median upper limit we expect to observe for a source lying in that particular 
declination. The current IceCube sensitivity is three times smaller than the 
previous IceCube 22-string sensitivity (Abbasi et al. 2009b). The upper 
limits set by IceCube are the best limits on neutrino astronomy so far. A 
source emitting an equivalent flux to the current sensitivity will be detected 
by the future IceCube in about 3-5 years, depending on the location in the 
sky. 



22 Strings 275.7 d 
40 strings disc. pot. 375.5 d 
40 strings sens. 375.5 d 
40 strings list 90% UL 
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Figure 5: Sensitivity and upper limits (90% C.L.) to a point-source E~ 2 flux 
as a function of declination. Blue solid line is the present median sensitivity based 
on this work. Green line is the previous sensitivity for 22-string configuration of 
IceCube analysis (Abbasi et al. 2009b, Abbasi et al. 2009c). Future expected 
sensitivity for IceCube, black solid line, and ANTARES (Coyle, 2010), dotted red 
line, are also shown. The squares shows the upper limits on the pre-selected 39 
direction of the sky for this work. 
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